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In this letter, it is reported that the first experimental results are associated with the GAM induced 
by energetic electrons (eEGAM) in HL-2A Ohmic plasma. The energetic-electrons are generated 
by parallel electric fields during magnetic reconnection associated with tearing mode (TM). The 
eEGAM localizes in the core plasma, i.e. in the vicinity of q=2 surface, and is very different from 
one excited by the drift-wave turbulence in the edge plasma. The analysis indicated that the eEGAM 
is provided with the magnetic components, whose intensities depend on the poloidal angles, and its 
mode numbers are |m/n|=2/0. Further, there exist intense nonlinear interactions among eEGAM, 
BAEs and strong tearing modes (TMs). These new findings shed light on the underlying physics 
mechanism for the excitation of the low frequency (LF) Alfvenic and acoustic fluctuations. 

PACS numbers: 52.35.Bj, 52.35. Mw, 52.35.Py, 52.35. Vd 
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Introduction-The very low-frequency (LP) Alfvenic 
and acoustic fluctuations, such as beta-induced Alfven 
eigenmode (BAE), and geodesic acoustic mode (GAAI), 
are presently of considerable interest in the present-day 
fusion and future burning plasmas [T], e.g. ITER. The 
low-frequency waves can significantly affect the plasma 
performance, and induce the particle losses and reduce 
the plasma self-heating. These LF instabilities can play 
an important role in turbulence and anomalous trans- 
port regulation, especially, while there is significant frac- 
tion of high energy particles in plasma |2] |3] . They can 
be used as energy channels to transfer the fusion-born- 
alpha-particle energy to the thermonuclear plasma, i.e. 
GAM/BAE channeling [3]. 

The GAM with toroidal mode number n=0 is an eigen- 
mode sustained by the coupling of radial electrostatic 
field and the poloidal variational density perturbations, 
and is usually taken to be electrostatic oscillation. The 
GAM is excited via modulation instability and pumped 
by the nonlinear interaction of drift wave turbulence [2] , 
and also driven by fast ions [S][B][7]. The GAM was in- 
vestigated both using gyro-kinetic simulations and ana- 
lytical methods in toroidal and slab geometry, and ob- 
served extensively in torus plasma [2][3]- Meanwhile, the 
BAE with n 7^ is also a low frequency mode with par- 
allel wave number fcii = (n — m/q)/RQ = , which is 
due to the plasma finite beta effect under the geodesic 
curvature, and usually believed to be electromagnetic os- 
cillation, and created by the coupling between the shear 
Alfven continuum with the poloidal mode number m and 
the sound continuum with the mode numbers m-1 and 
m+1, and driven by fast particles or large magnetic is- 



land. The BAEs were observed and investigated under 
different conditions in tokamak plasma "5^. 

It is worthwhile noting that the BAE and GAM have 
similar dispersion relations in the case of the long wave- 
length limit, i.e., the kinetic expression of the GAM dis- 
persion relation can degenerate with that of the LF shear 
Alfven accumulation point (BAE) [9] [10], which is use- 
ful for helping reciprocally identify the instabilities in 
the experiments. The most simple dispersion relations of 
BAE/GAM are given by 

L^BAE = OJGAM « (2T,/mO^/'(7/4 + TjT^y^yqRo 

Where q is safety factor, Rq is major radius, rrii is ion 
mass, and T^, T^ are ion and electron temperatures, re- 
spectively. The energetic-electrons and magnetic-island 
induced BAEs had been observed and investigated on 
HL-2A in the previous works [TT] [5]. In this letter, it is 
reported that the first experimental results are associated 
with the GAM induced by energetic-electrons (eEGAM) , 
and also present that there exists the intense nonlinear 
interactions among eEGAM, BAEs and strong TMs. 

Experimental conditions and mode characteris- 
t2cs-HL-2A is a medium-size tokamak with major/minor 
radius R/a — 1.65?7i/0.4?77,. The experiments discussed 
here were performed in deuterium plasmas with plasma 
current Ip ~ 150 - llQkA, toroidal field Bt ~ 1.32 - 
1.38r, and safety factor qa ~ 4.2 — 4.6 at the plasma 
edge. The line averaged density was detected by a hy- 
drogen cyanide interferometer. The polodial number m 
is measured using a set of seven Mirnov probes localized 
in the high field side (HFS) and eleven ones localized in 
the low field side (LPS). But the toroidal number n is 



measured using a set of ten Mirnov probes localized in 
the LFS of the vessel ^. Four CdTe scintillator de- 
tectors are placed outside the vacuum vessel in order 
to obtain information of the hard x-ray emission, and 
chordal distances of sight lines are rd=5, 9, 15 and 30 
cm, respectively. The range of the hard x-ray spectrum 
is E^ = 10 — 200keV divided into many energy bins by 
the PHA-software setting. 
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FIG. 1: Experimental parameters of the typical discharge 
with strong TM on HL-2A. Plasma current, Ip and density, 
He (a), magnetic probe signal (b), and corresponding spectro- 
gram (c), respectively. 

The eEGAM has been observed in the HL-2A Ohmic 
plasma for the first time, recently. This phenomenon is 
perfectly reproducible, and a typical discharge param- 
eters are shown in Fig[T] A coherent MHD fluctua- 
tion is visible around 17.5 kHz from 1250 ms to 2500 
ms. The toroidal mode number analysis indicates that 
this fluctuation does correspond to GAM due to n=0. 
In general terms, the magnetic component of GAM is 
two-order than the electric one, therefore it is very dif- 
ficult that it is observed in Ohmic plasma. However, 
the magnetic components of GAM had been observed 
in the same discharge. The analysis indicated that the 
poloidal number of GAM is m=2, and the fluctuation in- 
tensity depends on the poloidal angles. The phenomena 
can be interpreted by Zhou's theory [13] which suggests 
that the GAM has a magnetic component with m=2, 
which is created by the m=2 parallel return current, 
and the fluctuation intensity depends on the poloidal an- 
gles, i.e., Bg (X sin{26). The similar experimental results 
{Be oc sin(0)), which are associated with the density fluc- 
tuation induced by GAM, can be found in the previous 
document [T3]. The BAEs are also visible during strong 
TM activity with m/n=-2/-l in the same discharge. The 
characteristics of the BAEs were investigated in previ- 
ous works [8]. The mode numbers of the BAEs are 
m/n=2/l and -2/-1. There exists an island width thresh- 
old (~ 3.4cm) for the BAE excitation on HL-2A |H|. Note 
that the BAEs can not be completely explained by the 
present theory |15| . The magnetic fluctuation spectro- 



gram indicates that the GAM is always accompanied by 
strong TM and BAEs, and their frequencies comply with 

fa AM = fBAE2 — ItM, fcAM = fsAEl + JtM aS well 

as fa AM = {fBAE2 + /baei)/2. The GAM locahzes in 
the core plasma, i.e., in the vicinity of q=2 surface where 
the ion Landau damping 7^ is larger than the edge due to 
7i oc exp{—q^), and it is very different from one excited 
by the drift-wave turbulence in the edge plasma on HL- 
2 A [inj [17] . Such GAM is not observed in the absence of 
strong TM or BAEs. 
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FIG. 2: Enhancement of energetic electrons during magnetic 
reconnection at different CdTe channels on IIL-2A for shot 
#17455. Magnetic probe signal (a) and corresponding spec- 
trogram (e). Hard X-ray counts in arbitrary unit, (b)-(d) and 
(f)-(h). Left column, rd = 5cm; Right column, rd ~ 30cm. 
(b) and (f), E^ = 30- 40keV; (c) and (g); E^ ^ 40 ~ 50keV; 
(d) and (h), E^ — 50 — 60keV. Other energy bins do not been 
shown here. 

Relationship between energetic- electrons and 

EGAM-The existence of energetic-electrons during 
magnetic reconnection results in the excitation of GAM. 
Generation of energetic-electrons during magnetic re- 
connection has been the subject of a number of the- 
oretical and experimental investigations [T8] [ T9 ] [20] [ 2T] . 
The production rate depend critically on the ampli- 
tude of the electric field generated during reconnection. 
The electric field is JSy = {sBt/16rs)wmdwm/dt [IB] . 
where Wm is the width of the magnetic island, Wm = 
4:{BrrsRo/nsBty/^, and dwm/dt is the growth rate of 
magnetic island described by the tearing mode equation 
dwjn/dt — 1.2{r]/^Q)A'^ in the case of low beta. Here, 
A^ is the stability parameter, rj is the plasma resistivity, 
Ts is the radius of the magnetic surface, Br is the radial 
magnetic field perturbations, and s — {r/q)dq/dr is the 
magnetic shear. On the basis of experimental parame- 
ters, we can evaluate that electric fields are of the order 
of £^11 ~ 5V/m during the process of magnetic reconnec- 
tion on HL-2A. Analysis of HXR energy distribution has 



indicated that the energy of the energetic-electrons in 
flight is of the order of 20-200 keV. The time resolution 
of the PHA analysis did allow one to determine tempo- 
ral modifications of the spectrum. More details will be 
introduced in a separate paper. Fig[2] shows that the 
HXR fluxes with different energy bins increase with TM 
growing at t=1270 ms, and the eEGAM is also driven. 
Further, during strong TM, the energy distributions of 
energetic-electrons are all enhanced at different CdTe 
channels, shown in FigjSl and the non-Maxwell distri- 
bution beams exist in the core plasma, as a result, these 
energetic-electrons induce the excitation of eEGAM. 
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FIG. 3: Energy distributions of energetic electrons without 
and with BAEs and eEGAM at different CdTe channels on 
HL-2A for shot #17455. Black, blue and red lines are corre- 
sponding to t=1100-1110ms, 1300-1310ms and 1480-1490ms, 
respectively. 

Nonlinear interactions among BAEs, eEGAM 

and TMs-Thc nonlinear mode coupling can produce 
coherent mode structures which can provide overlap of 
wave-particle resonances in the minor radius, and trans- 
fer wave energy across different spatial scale. The role 
of nonlinear mode coupling is generally important in de- 
termining the mode excitation, saturation or damping. 
The nonlinear interaction also affects energetic particle 
redistribution/transport or plasma confinement. A novel 
result, which is nonlinear mode couplings among TM, 
BAEs and eEGAM, has been observed on HL-2A. For 
studying the nonlinear mode coupling, the squared bi- 
coherence [22| is given by 6^(/i,/2) =| BxYz{fi,h) \ 
I < I A(/i)y(/2) I' X I Z[h) f > with the Fourier 
bispectral BxYzifij2) =< X{h)Y{f2)Z*{h) >, h - 
/i ± /2 and < b\h,f2) < 1. Where X{f), y{f) 
and Z{f) are the Fourier transform of the time traces 
of a;(t), y{t) and z{t), respectively. The symbol <> de- 
notes the ensemble average over many realizations. It 
is convenient to represent the contribution of the non- 
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linear coupling from multiple modes to one mode with 
the summed squared bicoherence, which is defined as 
^bj,yz = ^}=f,±fJ>\hj2)/NU). Where N{f) is the 
numbers of realizations. To investigate the nonlinear 
mode coupling more precisely, the summed squared bi- 
coherence must be higher than the noise level, which has 
the value of (1 — bxyz/^)- Figffl shows the squared bi- 
coherence and summed squared bicoherence of a poloidal 
Mirnov signal. It is found that the nonlinear interaction 
between the fundamental ubae = 1 (or ubae = 1 ) BAE 
with Jbae and n = 1 TM with fxM at each different 
moment. The following matching conditions are satisfied 
among these modes, i.e. utm + nsAE 

Itm + Ibae = f BAE /GAM ^i' ^M and BAEs 

over, the n ± 1 BAEs interact with TM further and 

create a multitude of BAEs/GAMs, and I'bae/gam ^ 

ik-l)*fTM + fBAE, "^'bAE/GAM = C^- 1) *"TAf +"SA_E, 

where k is positive integer. If only considering the fun- 
damental frequencies of TM, BAEs and GAM, their non- 
linear interaction relations are /bab2 — Ibaei = '^Jtm, 
Ibae2 + Ibaei = '^Igam, and other expressions are 
fBAE = Igam ± Itm (i.e., GAM + TM ^ BAE) or 
Igam - Ibae ± /ta/ (i.e., BAE + TM ^ GAM). 
Need to point out here that the direction of wave en- 
ergy transfer is unknown, and it needs to be assessed. 
The auto-bicoherence of magnetic fluctuation is shown 
that I /2 1= fc * Jtm and \ f2±fi\^k* Jtm are strong 
coupling lines. It suggests that the nonlinear mode cou- 
pling process between TM and BAEs is similar with the 
coupling of drift- wave turbulence and GAM [2] [3] as well 
as TAE and EPM [23l. 
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FIG. 4; Squared bicoherence 6^ ^ ^ (fi, f2) and summed 
squared bicoherence E&^ a a (f) of Mirnov signal at t=1.4- 
1.5 s for shot #17455. Note that the resampling frequency is 
f=100 kHz, and nfft=512. 1/N denotes the noise level. 



Excitation Mechanism of eEGAM-The excita- 
tion mechanisms of eEGAM can be discussed, briefly. As 
we known, the energy transfer from energetic particle to 
the wave can be expressed by G oc u>df/dW + ndf /dP^p 



Here, df/dW and df/dP^ are, respectively, the parti- 
cle energy derivative of the distribution function of en- 
ergetic particle and the derivative for the toroidal mo- 
mentum. For the GAM, the driving free energy may 
directly come from the positive gradient of the distri- 
bution function df/dW and indirectly (because of n=0) 
from the radial derivative of distribution function by non- 
linear interactions. The magnetic island induces BAEs, 
then the eEGAM is excited via the nonlinear interac- 
tions among BAEs and strong TM which is a pump 
wave. While the matching conditions are satisfied, i.e. 
iL)i+L02 = 1^3 and ki + k2 = ks, the nonlinear interactions 
can occur between three waves in the plasma, and the 
coupled equations are bilinear and similar, and the cou- 
pling coefficients of the wave field amplitude determine 
the growth or damping of the waves [21] ■ According to 
Chen's theory, while the purely Alfvenic state described 
by the Walen relation {Su/va = ±B/Bo) is broken [25] . 
i.e. i5i?|| ^ or w^ 7^ ^n'^^Aj i^ ^'^^^ ^^^^ to significant 
perpendicular ponderomotive force and zonal flows. In 
our experiments, SE^^ 7^ is satisfied obviously owing 
to the magnetic reconnection of strong TMs. The slow- 
sound-wave density and potential perturbation, which 
are induced by parallel ponderomotive, will have radially 
varying (n = 0, m = ±1) poloidal structures. The par- 
allel ponderomotive force can couple with the compress- 
ible dynamics of slow-sound-waves, and the two high- 
frequency Alfvenic modes can generate a low-frequency 
acoustic mode by nonlinearity. If the perpendicular in- 
compressibility of shear Alfven wave is broken by the 
magnetic curvature, i.e. Su\ ^ 0, the AEs can nonlin- 
early generate (n = 0,?7i = ±1) radially local magnetic 
perturbations. Recently, the theory investigation sug- 
gests that the plasma compressibility has a significant 
effect on nonlinear mode coupling of AEs, and the cou- 
pling of AEs is a more effective energy transfer at a lower 
amplitude level due to the enhanced compressional per- 
turbations in the poloidal sidebands i26j. The GAM is 
toroidally symmetric mode unique to toroidal plasmas 
with mode structure that is nearly poloidally symmet- 
ric. The spatial overlap of the mode structures should 
be essential for the nonlinear mode coupling. Our ex- 
perimental results indicates the TM, BAEs and eEGAM 
all localize at the <\~2 surface, i.e., there exists an over- 
lap between the mode radial structures of eEGAM and 
BAEs. It means that the AEs can propagate poloidally 
into the region of the zonal fiows (ZFs) due to the zonal 
mode structure of GAM and the mode structure overlap, 
and can interacts with the GAM/ZFs, and the wave en- 
ergy can transfer between the GAM and AEs, and the 
intensity of the AEs can influence on the growth, satura- 
tion and damping of the GAM. 

Summary-The eEGAM has been recently observed 
for the first time on IIL-2A. The magnetic fluctuation 
spectrogram indicates that the eEGAM is always accom- 
panied by strong TM and BAEs. The analysis reveals 



that the eEGAM is provided with the magnetic compo- 
nents, whose intensities depend on the poloidal angles, 
and its mode numbers are | m/n |= 2/0. Further, a 
novel result, which is that there exist the cross-scale cou- 
plings among TM, BAEs and eEGAM, has been observed 
on HL-2A. The eEGAM is directly driven by energetic- 
electrons via the gradient of the velocity space or indi- 
rectly produced via the nonlinear mode coupling among 
BAEs and strong TM, but more theoretical works are 
needed because this phenomenon is a typical example 
with respect to multi-scale interactions. The eEGAM 
should have a significant effect on plasma transport in 
the vicinity of the magnetic island, and also have a pro- 
found regulatory effect on the turbulence around mag- 
netic island. The experimental results indicate that the 
couplings possibly induce the energy transfer among TM, 
BAEs and eEGAM, and it is possible to be one of mech- 
anisms of the energy cascade in Alfven turbulences, and 
the BAE/GAM may be an energy channeling between 
different scales, such as macro-, meso- and micro-scale. 
The new findings give a deep insight into the under- 
lying physics mechanism for the excitation of the LF 
Alfvenic/acoustic fluctuation and ZFs. 
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